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This article presents the design control of a tracking problem for a non- 
minimum phase bilinear control system containing disturbance. The bilinear 
control system is assumed to have a relative degree one and non-minimum 
phase, which means it has unstable internal dynamics. The disturbance exists 
only in state variables corresponding to the control function in external 
dynamics. The control design was carried out using the backstepping 
method, which was applied to the normal form of the bilinear control 
system. Internal dynamics will be stabilised using virtual control to 
overcome unstable internal dynamics. The last step will stabilise the external 
dynamics and disturbance using the original control function. The simulation 
results show that the proposed control method can rapidly drive the output to 
the given trajectory. Control performance varies depending on the control 
parameter setting. The higher the control parameter, the better the control 
performance, evaluated using integral absolute error. 
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1. INTRODUCTION 

The dynamical system is widely used in modelling and analysing natural phenomena, both physical 
and social behaviour. The accuracy in the analysis process is highly dependent on the modelling process, 
including the selection of the type of model used [1], [2]. Models with nonlinear equations are the best form 
of modelling, but nonlinear systems are complicated to analyse, so approximations are often made to 
interpret them. The bilinear control system is one of the classes and is the best approximation of the nonlinear 
control system [3], [4]. The approximation of a nonlinear control system using a bilinear control system can 
be made using the Carleman transformation [3] or the Jacobian transformation [5]. The advantages of a 
bilinear control system in approaching a nonlinear control system compared to a linear control system in 
terms of performance, optimisation and modelling processes [6]-[8]. Intensive studies on the use of bilinear 
control systems in several fields in economics and chemistry can be referred to [1], [9], [10]. 

Some examples of the use of other bilinear control systems can be seen in the boost converter 
problem in [11], [12], on controlling the spread of disease in [13], [14], on directing ship motion in [15], on 
regulating hot water storage in [16], on the diesel engine fuel collector in [17]. The application of control to 
the Hilbert chamber using finite feedback has also been carried out [18]. The model's accuracy is also 
determined from the modelling process carried out. According to Amato [19] the dynamical system contains 
an uncertainty factor due to simplifying the modelling and parameter estimation processes. The first 
uncertainty factor occurs in the modelling process because the phenomenon's complexity cannot be fully 


Journal homepage: http://ijeecs.iaescore.com 


1316 O ISSN: 2502-4752 


expressed in a mathematical model and natural phenomena contain a constant disturbance function that 
cannot always be modelled [2]. 

Inaccurate estimation of model parameters also results in an uncertainty system. The dynamical 
system, which represents the phenomenon, will contain parameters whose values cannot always be 
determined precisely or whose values change at certain time intervals. The uncertainty factor in a 
mathematical model can be expressed in the form of a function that represents a parameter whose value is 
constant but unknown, a time-varying parameter, or a function that defines an external disturbance [21]. One 
of the problems arising from a disturbance function is not achieving asymptotic stability from the system 
output [22]. 

One of the methods used in the nonlinear control system is the backstepping method. This method 
was developed in the 1990s based on the stability of Lyapunov which the process is carried out iteratively 
[23]. Using the backstepping method can be referred to in [24]-[29]. Generally, the backstepping method is 
applied to systems with a strict feedback form [30]—[32], which can be obtained using the transformation in 
[33]. This transformation is a nonlinear transformation that can make the control design more complex. 

This article presents the control design of a non-minimum phase bilinear control system containing a 
disturbance function. The main difference with previous studies in [12]-{14], [16], [18] is in terms of the 
presence of a disturbance function in the system and the type of system and in [15], [34]-[36] is in terms of 
the use of the control design method used. The method used in this article is the backstepping method, which 
develops the previous results in [37]. In [37], it is assumed that the system can be linearised exactly so that it 
has no internal dynamics, and the control design is only a matter of stabilisation. Compared with the work in 
[27]-[29], the system used is not in a non-minimum phase. Even though the system contains uncertainty, the 
absence of internal dynamics makes it possible to design controls directly on the system. We applies 
backstepping to the normal form obtained through the input-output linearisation transformation using the 
procedure in [38], [39]. Since there are internal dynamics, they will be stabilised first using virtual controls. 
The whole system and the disturbance function will be stabilised in the last iterative process using the actual 
control variables. The novelty aspect of this article includes theories and methods in solving a non-minimum 
phase system containing a disturbance function. To solve the non-minimum phase system, it is no longer 
necessary to redefine the system output to become a minimum phase. In addition, the use of backstepping 
ensures that the internal dynamics can be stabilized to zero or bounded using virtual control. 

This article is structured as: After this sections, a research method explains the input and output 
linearisation transformations applied to a bilinear control system. It also describes the basic backstepping 
procedure to stabilise the bilinear control system and the method for calculating control performance using 
integral absolute error (IAE). In the results section, the control design for the non-minimum-phase bilinear 
control system containing disturbance is described and analysed for its stability. Some examples are provided 
to simulate control design and to see the performance of the proposed method using several different control 
parameters. 


2. RESEARCH METHOD 
Given a bilinear control system with a single input and output: 


es 
y(t) 
where A,B E€ R2*?, x(t) € R? is the state variable, u(t) € R is the control variable, and y(t) is the system 
output. It is assumed that the system output can be expressed as a linear combination of state variables, i.e., 
h(x(t)) = hx(t) for h” € R?. 

Linearisation of the bilinear control system is done based on input-output feedback using the Lie 
derivative defined by Lph(x(t)) = Vh(x(t))f (x(t). The relative degree of a bilinear control system in (1) 
is the natural number p € N that satisfies Li*h(x(t)) = hA'1Bx(t) = 0 for each i = 1,2,...,9 — 1 and 
Lgl *h(x(t)) = hA?1Bx(t) #0. If p<n , then a bilinear control system has internal dynamics 
independent of the control variables and if p = n is satisfied, then a bilinear control system in (1) is a system 
that can be linearly exactly. If the system's internal dynamics are unstable, then the system is called a non- 
minimum phase system. 

Assuming the bilinear control system in (1) has a relative degree two which satisfies Lgh(x(t)) =0. 
The input and output linearization transformation (IOFL) transformation are given by (2). 


Ax(t) + u(t)Bx(t) 
h(x(t)) 


(1) 
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h 
E a 2) 


The derivative of (2) with respect to t and using its inverse gives: 


ee = z(t) 


z(t) = v(t) @) 


where v(t) = hA?T~1z(t) + u(t)hABT~1z(t) is the new control function that applies to the linear control 
system in (3). Choose Lyapunov function V(t) = zz? + : (Z> + r1z1)?. The linear control system in (3) can 
be stabilised using the backstepping method with the control function v(t) = (r? — 1) z — (r +r) x for 
1 2 
Ty. E Rt. 
It is assumed that the bilinear control system in (1) has a relative degree one to satisfy L,h(x(t)) + 
0. Defined coordinate transformation using input and output linearization transformation (IOFL) [38], [39]. 


h(x) 
P(x) 


Where (x) is a function that satisfies L(x) = 0 with g(x) = Bx(t). The derivative of (4) with respect to 
t and using the inverse transformation of (4) produces: 


z=T(x) = | 4) 


o = v(t) 6) 


Z(t) = Leb(x) = c12, (t) + c222(t) 


where v(t) = RAT~1z(t) + u(t)hBT~1z(t) is the new control function that applies to the linear control 
system in (5). In contrast to the system that can be exactly linearised in (3), the control value does not affect 
internal dynamics. If the controlled output value reaches the path, then the dynamic zero on the internal 
dynamic is Z,(t) = c,z(t). If cp E R*, then dynamic zero is unstable so that the system in (1) is the non- 
minimum phase. 

The control function parameter {r,,7,} in the backstepping method determines the control 
performance to bring the system output to follow the given path. Performance values are evaluated using the 
integral absolute error (LAE), which is defined by (6). 


IAE = f ly(t) — ya(t)| dt (6) 


The integral in (6) needs the numerical solution of the system of (1) so that the value of (6) is evaluated using 
Simpson's method given by; 


h = n 
Sy = (fo oe ees E eee 232 9 fas] 


tr—to 


where n is the number of partitions used, h = is the width of the partition, t; = tọ + ih is the partition 


point for the t domain and f; = |y(t;) — yq(t,)| for i = 0,1,2,...,n. 


3. RESULTS AND DISCUSSION 
This section presents a control design for tracking problems in a bilinear control system with a 
disturbance function using the backstepping method for a non-minimum phase system. Simulation examples 
are also given to see the performance of the control in bringing the output along the given trajectory. 
Consider a planar bilinear control system with a disturbance function. 


o 
y(t) 


Ax(t) + u(t)Bx(t) + Nw(t) 


h(x(t)) 0) 
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where A, B € R2*?, x(t) € R? is the state variable, u(t) € R is the control variable, and y(t) = h(x(t)) is 
the system output assumed to be expressed as a linear combination of the state variables, so that y(t) = 
hx(t) and w(t) € R are disturbance functions. 

- Assumption 1. The tracking path y,(t) and its derivatives are smooth and bounded. 

- Assumption 2. The disturbance function w(t) is unknown but bounded. 


3.1. Control design for the system with exact linearisation 
Theorem 1, given a bilinear system with disturbance in (7). If the system has a relative degree of 
v(t)—-hA? x(t) 


FABRE. 2 where v(t) is given by: 


two, the control function that stabilises the system globally is u(t) = 


v(t) = Jal) — (1 — rf es (t) — (rı + rw)w(t) — sign(w) max(d) 


Proof, it is known that the system has a relative degree of two to be linearised exactly. Using the 
transformation of (2), we get (8): 


Z(t) = v(t)+d(t) 
defined the difference between the output and the given path (9), 
fa) = a-n ies 
e(t) = 22(t)— ya(t) 
the derivative of (9) with respect to t and by using (8) gives (10), 
ee = e(t) (10) 
êt) = v(t)+d(t)— Yalt) 


where d(t) = hNw(t) is the disturbance function and v(t) is the new control variable used in the linearised 
system of (10), and its relationship to the control variable u(t) is: 


v(t) = hA?x(t) + u(t)hABx(t) (11) 


The first step is to stabilise the state variable e,(t) using e,(t) as a virtual control. Define the 
Lyapunov function V; (t) = ~e?(t). The derivative of V; (t) with respect to t gives: 


V(t) = e, (t)e2(t) 


For the state variable e,(t) to be asymptotically stable, V, (t) < 0 must be valid for every t. Select virtual 
control e,(t) = —r,e,(t) with r, E€ R*. Define a new state variable (12): 


w(t) = rie, (t) + e2(t) (12) 
based on (12), we get e(t) = w(t) — r,e,(t), and by substituting in (10), we get (13). 
ė (t) 
w(t) 
The final step is to stabilise the entire system in (13) using the variable control v(t). Define the 
Lyapunov function V(t) = 5 [e?(t) + w? (t)]. The derivative of V, (t) with respect to t will results, 


w(t) — r,e, (t) 
rı (wlt) — rie lt)) + v(t) + d(t) — Yalt) 


(13) 


VE) = -re (t) + w(t)[qw(t) + A - rPe) + vE) = Pat] + wO) (14) 


It is assumed that there are r, € Rt and ® € R™ such that V(t) = —r,e?(t) — r„w?(t) — ®. Using this 
assumption and using (14), it obtained in (15): 
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r,w(t) + (1 —7Z)e,(t) + v(t) — Yalt) = —n, w(t) + 0. (15) 
Substitute (15) into the Lyapunov function in (14) to get (16). 
V2(t) = —r,e? (t) — ryw? (t) + w(t) (8 + d(t)) (16) 


To satisfy V,(t) < 0, select 0 = —sign(w) max(d). Substitute @ = —sign(w) max(d) into (16) so that the 
final form of the derivative of the Lyapunov function is obtained that is (17). 


V2 (t) = —rye7 (t) — rww? (t) + |w|(sign(w)d(t) — max(d)) (17) 


Because max(d) > d(t) then sign(w)d(t) — max (d) < 0 and V,(t) < 0 apply to every t > 0. 
Furthermore, the variable control v(t) is obtained from (15) and using the relation between the state variables 
{e,(t) , e2(t) , w(t)} in (12), the following control function is obtained (18). 


v(t) = Yalt) — (1 — rf )er(t) — (rı + rw)w(t) — sign(w) max(d) (18) 


Using the Lyapunov function V(e,,w) and its derivatives in (17), we found that the value {e,(t),w(t)} is 
bounded. Use Assumption 1, and because the value of yg(t) is bounded, based on (9), the value of z,(t) is 
also bounded. Because {e,(t),w(t)} are bounded, from (12), the value of e, (t) is bounded. As a result, from 
(9), we obtained 2,(t) is also bounded. Based on the LaSalle- Yoshizawa theorem, we get lim (e(t)) = 0, so 


lim (y(t) — ya (t)) = 0 applies. 


3.2. Control design for non-minimum phase system 
Theorem 2. Consider a bilinear control system in (7) with relative degree p = 1 and satisfies 
Assumption 1 and 2. Let [ky k2] = cAT~1 with T = (h,c) and c is an internal dynamic equation matrix 


that satisfies cBx(t) = 0. The control function that takes the output to the path yg(t) is u(t) = a 


where the value of v(t) is given by: 


M'Tx(t) 
ky 


v(t) = Ya CE) + (ke + Ta + rw) yal) — ES + E MECE) — SI sign (w) max(d) 


ka 


Whee Nd eae es bE) 


and the function @(t) are bounded solutions of the first-order 


differential: 


p(t) = k(t) = kia (t) =0 


Proof, it is known that the bilinear control system of (7) has a relative degree p = 1 so that it has 
internal dynamic and internal dynamic stabilisation will be performed first. Using the transformation in (4), 
we get (19). 


ae = v(t)+d(t) (19) 
Z(t) = kız, (t) + k222(t) 
Define the difference between the output and the given path, 
e(t) = z4) - yalt) 

20 
oE O e 


where (t) is a function to be defined later. The derivative of (20) with respect to t and using (19) produces 
(21). 


Bee = k,(e,(t) + ya(t)) + ko(e2(t) + (0) - 6) ei 


élt) = v(t)+d(t) -— yalt) 


The first step is to stabilise the state variable e(t) using e(t) as a virtual control. Define the 
Lyapunov function V(e,) = ~e3(t). The derivative of V (e3) with respect to t produce: 
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V (e2) = e2(t)(ky(e,(t) + ya(t)) + kp (e2(t) + o(t)) - 940) 


the state variable e,(t) is asymptotically stable if it satisfies V(e,) < 0 for every t > 0. Choose virtual 
control k,(e,(t) + Yq (t)) +k, (e2(t) + p) — P(t) = —re,(t) with rı € R*. Define a new state 
variable (22). 


w(t) = kyex(t) + (kz + r2)e2(t) + kiyalt) + k24) — b(t) (22) 
Substitute new state variable in (22) and its derivative into (21) to produce (23). 


K (t) w(t) — me2(t) s 
w(t) = ky(v(t) + d(t) — ya(t)) + (kz + )(W(E) — rzez (t)) + kida lt) + kz$ (t) — pC) 


(23) 


The final step is stabilising the entire system in (23) using the variable control v(t). Define the 
Lyapunov function V(e2,w) = = (ef (t) + w?(t)). The derivative of V(e2,w) with respect to t will results 
(24). 


Ý (e2,w) = =rze3 (t) + w(t) (e2 (t) + w(t) (24) 


Substituting (23) into (24) and assuming there are n, € Rt and ® € R™ such that V(e2,w) = —r,e3(t) — 
r,W2(t) — ® to get (25). 


e(t) + ky (v(t) — Yalt)) + (kz +) (W(t) — re2(t)) + ky Walt) + koh (t) — H(t) 
= —tyw(t) + 0 (25) 


Using (25), (24) became (26): 
V(e2,w) = —r,e3(t) — rw? (t) + w(t)(k,d(t) + 8) (26) 


To satisfy V(e,w) < 0, choose the value @ = —sign(w)|k,|max(d). Substitute this value into (26), and we 
obtain the final form of the derivative of the Lyapunov function, that is: 


V(e2,w) = =r2e3 (t) — row? (t) + |w(t)k,| (sign(k,)sign(w)d(t) — max(d)) (27) 


Because max(d) > d(t), then sign(k,)sign(w)d — max(d) < 0 and V(e2,w) < 0 apply to every t > 0. 
The control function v(t) is obtained from (25). 


v(t) = -> 


[(1 — r2(k2 + r2))e2(t) + (k2 + r2 + ty)w(t) — p(t) + kz$ (Ct) + sign(w)|k; | max(d)| (28) 


Using the value of {e4 (t), e2(t)} in (20) and w(t) in (22), the control function v(t) in (28) can be expressed 
by (29): 


M'Tx(t) 
ky 


v(t) = Ya CE) + (ke + Ta + tw) yal) — EO + E MGC) — SI sign (w) max(d) (29) 


ka 
: 1. ky (kz + r2 + Ty) 
Wak (este, Week HOF 

Note that the control function in (29) requires a value of @(t). Consider the w(t) function in (22); 
w(t) = k,e,(t) + (kz + r2)ez2(t) + kyyg(t) + k2b(t) — (t). The state variables e(t) and w(t) have 
been controlled to zero at each stage. Furthermore, it is expected that e,(t) ~ 0 to obtain z,(t) = y(t). 
Assume that for a value t > T that satisfies e,(t),w(t) ~ 0, from (22), we will get k,e,(t) = $(t) — 
kzp (t) — kiya (t). Next, to obtain e(t) = 0, the value of @(t) is obtained from an ordinary first-order 
differential equation. 


A(t) — kzt) — kiya lt) = 0 (30) 
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The solution of (30) is a path for internal dynamics, and it depends on the initial value of (tọ) and 
the path yg(t). In order to make bounded the internal dynamic, it is necessary to determine the initial value 
of $ (to) so that the solution of (30) is bounded for every t > 0. 


3.3. Simulations 

This section gives some examples to illustrate the implementation of the controls described in the 
previous section. Example 1. Consider the dynamics of the paper cutting machine [34], which is expressed by 
(31): 


x(t) = —0.046x,(t) — 0.027x,(t)u(t) + 0.978d j 
P = —0.7632x,(t) + 3.197x2(t) a) 
Suppose v(t) = —0.046x,(t) — 0.027x,(t)u(t) then (31) can be written as (32). 
ha = v(t) +0.978d 39 
p(t) = —0.7632x,(t) + 3.197x»(t) 2) 


As shown in (32) has formed a normal system with a control function v(t) with an internal dynamic of order 
one. If x, = 0, then zero dynamics is x(t) = 3.197x,(t), which is unstable and consequently is a non- 
minimum phase system with relative degrees p = 1. Simulations were carried out using the data in Table 1. 
Using Table 1, the control function in (18), which takes the output to the origin, is v(t) = 
—1.3103(17.5910x, — 11.1970w — 0.7632 sign(w)). Figure 1 shows the simulation results for the 
stabilisation problem. Figure 1(a) shows that the proposed control function with several parameter values 
variations successfully stabilises the system output to the origin. The system output initially moves from the 
initial value away from the origin to a certain maximum point and then towards the origin. The convergence 
behaviour of the system output is strongly influenced by the value of the backstepping control parameter. The 
maximum point achieved for a small combination of control parameters is smaller, but the system output 
takes longer to converge to the origin. As for larger combinations, it takes a relatively short time to converge, 
but the system will move to a larger maximum point. Figure 1(b) shows the internal dynamics for each 
backstepping control parameter. In contrast to the system output, the internal dynamics controlled with large 
parameters move from the initial value to the origin directly. Meanwhile, the internal dynamics, which are 
controlled with small parameters, allow the internal dynamics to move away from the origin until a certain 
maximum point before moving towards the origin. 


Table 1. Simulation parameters 


Parameter Value Description 
T2 {3,7,11} Backstepping parameter 
Ty {5,9,15} Backstepping parameter 
d(t) rand() X sin(t) Disturbance 
kı —0.7632 Coefficient variable internal dynamic 
k2 3.197 Coefficient variable external dynamic 


Next up for the tracking problem, we use the path yg(t) = cos(t). The path for internal dynamics 
obtained from the equation @(t) is obtained from the solution of the following ordinary differential equation. 


49) _ 3.1976(t) + 0.7632 cos(t) = 0 (33) 


dt 
If the initial value is (to) = @o , then the solution of (33) is (34). 


p(t) = 0.217 cos(t) — 0.068 sin(t) + exp(3.197) (po — 0.217) (34) 


In order for (34) to be bounded, the term exp(3.197t) needs to be eliminated. Selecting the initial value of 
o = 0.217 will produce a bounded solution of (34) as well as a particular solution, namely (t) = 
0.217 cos(t) — 0.068sin(t). The simulation result for the tracking problem with the path yg(t) = cos(t) is 
shown in Figure 2. Figure 2(a) compares the system's output dynamics with the trajectory. The system output 
has converged to the trajectory at t ~ 1. At t € [1,10], the system's output and the trajectory remain 
attached, and there is no dynamic spike even though the system is disturbed. Figure 2(b) shows the internal 
dynamics. Internal dynamics move from the initial point following a path @(t) whose value is bounded at 
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[—0.5,0.5]. This indicates that the backstepping control is successful in bringing the system output to the 
trajectory while keeping the internal dynamics to a minimum. 


(0), (r2 = 3,r» = 5} 
-= -x1 (t), {re = T, Tw = 9} 
TRE z(t), {re = to = 15} 


1.2 T T T 
x(t), {r2 > 3, Tw = 5} 
i - = =29(¢), {r2 = T, ru = 9} l 
RN O peee x(t), {r2 = 11, ru = 15} 
0.8} 3) 
ia 
ia 
PE aE 
5 
0.47 
0.27 
ob 
-0.2 " 
0 0.5 1 5 2 2.5 3 


Figure 1. Comparison of state variable dynamics on paper cutting machines with varying control parameter 
values with (a) system output and (b) internal dynamic 


Example 2 given a bilinear control system with a disturbance function x(t) = Ax(t) + u(t)Bx(t) + 
Nw(t) with the system output y(t) = x2(t) and the parameter values are given in [40]. 


e a aal®=lo e G5) 


The disturbance function is generated using the equation w(t) = rand( ) X sin(t). The purpose of the 
control is to drive the system output to the origin yg(t) = 0. For the control parameters, the values of {r, = 
1,%, = 3} and {r, = 5,7, = 7} are used. Comparison of computational results between the proposed method 
and the robust control method Ho [13], [40] is shown in Figure 3. 

Figure 3 shows that both control methods can bring the output to the origin relatively quickly. 
Compared with the robust Hoo control, the results depend on the combination of {1,,7%,} values. In the 
backstepping method, the control performance shows significantly different results for the two variations of 
the given {r}, ny} value. For the {r, = 1,7, = 3} combination, the backstepping method produces worse 
performance than the Hoo robust control method. The system output controlled using backstepping takes 
longer to converge to the origin. Still, the robust Hoo control method produces a higher minimum output 
value at the beginning of the simulation. By changing the value of the control parameter {r,,7,} to be larger, 
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the convergence of the backstepping method increases rapidly. The {r, = 5,7, = 7} combination can bring 
the system output with a minimal minimum value to the origin in time t ~ 0.5, which is 700% faster than the 
{ri = 1,% = 3} variation, which takes time t ~ 3.7. 


| =a 
- - -ya(t) = cos(t) 


Figure 2. Comparison of state variable behaviour in the tracking problem with (a) system output versus the 
path and (b) internal dynamic 


For tracking problems, select the path ya (t) = cos(t). The implementation of the control function 
in (29) using several variations {r;, rw} is shown in Figure 4. The calculation of the error value between the 
system output and the path using IAE for each {r}, ny} combination is shown in Table 2. 

Figure 4 shows, with all given combinations value of {r}, y}, the system output can follow the path 
with identical behaviour but at different speeds. At the start of time, the system output moves away from and 
against the given path. The system output moves back to the given path at a maximum point, shown at time 
t ~ 0.7. The smaller the combination of {r;, ry} values, the higher the system output value and the longer it 
takes to converge to the path. These results provide a different IAE value for each combination {r,,7,,} with 
the characteristic that the lower the combination value, the higher the IAE value for the combination. Table 2 
shows the IAE values of each combination shown in Figure 4. The combination {r, = 3,n, = 3}, the 
smallest combination, gives the largest IAE value, while the combination {r; = 12,7, = 9}, the largest 
combination, gives the smallest IAE value. 
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Pre - - -Backstepping, {r;=1, r,=9} 4 


Robust 


Backstepping, {r1=3, ru=13} 


4 5 6 7 8 9 10 


t(s) 


Figure 3. Comparison of the system output stabilised to the origin using the backstepping method with two 
variations of {7,,7%,} and the robust Hoo control method 


-—--y(t),{r1=5,rw=7} 
= u(t) {ri=7ru=9} 
— y(t), {r1=12,r,=9} 
—-—-— yalt) = cos(t) 


4 5 6 7 8 9 10 


Figure 4. The effect of variations in the value of {r}, ny} on control performance 


Table 2. IAE values based on the combination {r,,7,,} shown in Figure 4 


{rs fw} IAE ToC (seconds) 
M =3,% =3} 7.63929 0.58891 
(1, =5,%y =7} 1.96321 0.54580 
{ri =7,ůy = 9} 1.04368 0.54813 
{r, = 12,7, = 9} 0.60237 0.54495 


Table 2 does not adequately explain the change in IAE value due to changes in the value of the 
combination of {7,,7%,}. Next, the IAE calculation is performed using a combination of partitioned (r4, ny) = 
[1,25] x [1,25] with Ar = 0.5, and the results are shown in Figure 5. The highest IAE value occurs for the 
combination of {r,; = 1,%, = 1}, and the lowest IAE value occurs for the combination of {r, = 25,n, = 25}. 

Statistically, the correlation level can be calculated from the combination {r;, rẹ} to the IAE value. 
Using Spearman's non-parametric correlation test, r, to IAE is —0.68043 while n, to IAE is —0.68240 with 
1% significance level. This means that {r;, ny} has enough correlation to lower the IAE value. Furthermore, 
if multiple regression is performed to obtain the IAE function for the {r,,7%,} combination, the equation 
IAE = —0.106467, — 0.106577, + 3.64456 is obtained. Although the weights of the two coefficients 
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{r1,%y}are different, the difference is very small. We say that both {r,,%,} give almost equal weight to 
changes in IAE values. 


Figure 5. IAE value for all combinations of (r4, ry) = [1,25] x [1,25] 


4. CONCLUSION 

This article discussed the control design for tracking problems using the backstepping method for a 
non-minimum phase bilinear control system containing disturbance. The disturbance function is assumed to 
exist only in dynamic externals. The bilinear control system is converted into a normal form using a 
linearisation transformation based on input and output. The control function is designed using the 
backstepping method from the normal form. For systems that have internal dynamics, the internal dynamics 
are stabilised first. For the tracking problem, a path has been defined for the internal dynamics so that the 
system output follows the path, and the internal dynamics remains bounded. Based on the simulation, the 
proposed control function successfully brings the output to follow the trajectory for stabilisation problems 
and tracking problems. The selection of the right parameters will affect the performance of the control and, 
for a certain value, results in better performance than the Robust method on stabilisation problems. The 
combination of control parameters is almost symmetrical and inversely proportional to the resulting IAE 
value. 
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